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ABSTRACT 

We examine unresolved nuclear X-ray sources in 57 brightest cluster galaxies to study the 
relationship between nuclear X-ray emission and accretion onto supermassive black holes. 
The majority of the clusters in our sample have prominent X-ray cavities embedded in the 
surrounding hot atmospheres, which we use to estimate mean jet power and average accretion 
rate onto the supermassive black holes over the past several hundred Myr. We find that roughly 
half of the sample have detectable nuclear X-ray emission. The nuclear X-ray luminosity is 
correlated with average accretion rate determined using X-ray cavities, which is consistent 
with the hypothesis that nuclear X-ray emission traces ongoing accretion. The results imply 
that jets in systems that have experienced recent AGN outbursts, in the last ~ 10 7 yr, are 'on' 
at least half of the time. Nuclear X-ray sources become more luminous with respect to the 
mechanical jet power as the mean accretion rate rises. We show that nuclear radiation exceeds 
the jet power when the mean accretion rate rises above a few percent of the Eddington rate, 
or a power output of ~ 10 45 erg s _1 , where the AGN apparently transitions to a quasar. The 
nuclear X-ray emission from three objects (A2052, Hydra A, M84) varies by factors of 2 — 10 
on timescales of 6 months to 10 years. If variability at this level is a common phenomenon, it 
can account for much of the scatter in the relationship between mean accretion rate and nu- 
clear X-ray luminosity. We find no significant change in the spectral energy distribution as a 
function of luminosity in the variable objects. The relationship between accretion and nuclear 
X-ray luminosity is consistent with emission from either a jet, an ADAF, or a combination 
of the two, although other origins are possible. We also consider the longstanding problem 
of whether jets are powered by the accretion of cold circumnuclear gas or nearly spherical 
inflows of hot keV gas. For a subset of 13 nearby systems in our sample, we re-examine the 
relationship between the jet power and the Bondi accretion rate. The results indicate weaker 
evidence for a trend between Bondi accretion and jet power, primarily due to the uncertainty 
in the cavity volumes. We suggest that cold gas fuelling could be a likely source of accre- 
tion power in these objects however we cannot rule out Bondi accretion, which could play a 
significant role in low power jets. 

Key words: X-rays: galaxies: clusters — galaxies: active — galaxies :jets — accretion, accre- 
tion discs 



1 INTRODUCTION 

Energetic feedback from supermassive black holes (SMBHs) plays 
an important role in the formation and evolution of galaxies. 
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With this realisation, the exploration and understanding of ac- 
tive galactic nuclei has taken a new emphasis. Key lines of ev- 
idence include the relationship between nuclear b lack hole mass 
and t h e mass of the host galaxy (M - o relation; [M agorrian et al.1 
1 19981 : iKormendv & Gebhardtl 1200 ll : iMerritt & Ferraresd 1200 lE 
which was likely imprinted th r ough the quasar era ('quasar mode 
feedback' : Silk&Rees 1998; Haehnelt Natarajan & Rees 1998; 
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Di Matteo. Springel & Hernquistll2005h . and the prevalence of ra- 
dio bubbles in the X-ray atmospheres of gi ant elliptical and bright- 
est cluster galaxies ('radio mode feedback' ; IMcNamara et al.l2 000 ; 
Churazov et all l2000l : iFabian et alJ l2000l : IMcNamara & Nulsenl 
2007h . These forms of feedback are ultimately powered by binding 



energy released by accretion onto massive black holes. However, 
energy is released primarily in the form of radiation from quasars, 
while radio jets release their energy in a mechanical form. 

The reason why black holes release their binding en- 
ergy in different forms is poorly understood. Clues have 
come from the so-called f undam e ntal plane of black holes 
(IMerloni. Heinz & di Matted 120031 : iFalcke. Kording & Markoffl 
This relationship between the radio and X-ray power emerg- 
ing from the vicinity of the black hole and the black hole mass ex- 
tends over nine orders of magnitude from stellar mass black holes 
to supermassive black holes. The continuity of this relationship 
indicates that the emergent properties of black holes are funda- 
mentally the same regardless of their mass. Furthermore, accreting 
black hole binaries undergo changes in emission states that seem 
to correlat e with the accretio n rate normalized to the black hole 
mass (eg. iFender et al.1 Il999|: iMaccarone. Gallo & Fended 120031 : 
iRemillard & McClintockll2006) . 

Changes in the specific accretion rate apparently lead to struc- 
tural changes in the accretion disk which govern the release of 
binding energy in the form of a jetted outflow or radiation. Dur- 
ing periods of high accretion rate a geometrically thin and opti- 
cally thick disk fo rms that dissipates its ene r gy primarily in the 
form of radiation dShakura & Sunv aev 1973; iNovikov & Thornd 
1 19731 : iFrank. King & Raindl2002h . During periods of more mod- 
est accretion, a hot, radiatively inefficient accretion flow (RIAF) or 
advection-dominated accretion flow (ADAF) forms t hat releases its 
energ y in the form of a jetted ou t flow or wind (eg. iNarayan & Yil 
1994; lAbramowicz et all 1 19951 : iNaravan & McClintockl 12008). 



ADAFs are geometrically thick, optically thin disks where the 
ion temperature substantially exceeds the electron temperature. In 
these disks, the inflow timescale is much shorter than the cooling 
timescale. Thus the accretion energy cannot be radiated and is ei- 
ther advecte d inward with the flow or is re l eased in a wind or radio 
jet (ADIOS; Blandford & Begelman 1999. 120041) . In the context of 
SMBHs, a radiatively efficient disk is formed when the accretion 
rate approaches the Eddington value giving rise to a Seyfert nu- 
cleus or quasar. When the accretion rate falls below a few percent 
of the Eddingto n rate, the nucleus be comes faint and a radio galaxy 
is formed (eg. Churaz ov et al 1l2005h . 

Galaxy formation models incorporating AGN feedback dis- 
tinguish between radiatively-dominated quasar feedback at early 
times and a mechanically-dominated radio mode at late times 
(eg. ISpringeL Di Matteo & Hernauisil 120051: ICroton et al I Exffi 
ISiiack i & Springel 2 0061 : iHopkins et al.| [2006: Bo wer et al.ll2006h . 
Quasar feedback operates through intense radiation that is expected 
to couple to the gas and strong winds which drive gas from the host 
galaxy, quenching star formation and regulating the growth of the 
SMBH. This will eventually starve the SMBH of fuel and, as the 
accretion rate drops, a transition to mechanically-dominated radio 
mode feedback is expected. The SMBH launches jets which reg- 
ulate radiative cooling in the surrounding hot atmosphere and the 
growth of the most massive galaxies. The AGN activity is closely 
correlated with the prope rties of the host halo indicating that they 
form a feedback loo p (eg. lBirzan et al. 2004; Dunn & Fa bian 2004; 
iRaffertv et ap2006h . 

Despite evidence for rapid accretion onto their SMBHs, gi- 
ant elliptical and brightest cluster galaxies in the nearby Universe 



rarely harbour quasars. The exceptions include H1821+643, 3C 186 
and IRAS09104+4 1 09 ([Crawford et al.lll999l: iBelsole et alJl200l 
iRussell et all l2010l : ISiemiginowska et all I2OI0L Cavagnolo et al. 
submitted). Instead most harbour low luminosity AGN ie. ADAFs 
(IFabian & Reeslll995l:ISambruna et al.l2000l : lDi Matteo et al.l2000l . 
Hlav acek-Larrondo & Fabianll201ll) . These low luminosities imply 
that their host galaxies harbour massive black holes exceeding ~ 
1O 9 M0. Nuclear emission likely signals ongoing accretion. How- 
ever, any dependence of the amplitude and form of the emission 
emerging from the nucleus on the accretion rate, as found in X-ray 
binaries, requires an independent means of estimating the accretion 
rate itself. Therefore it is difficult to test nuclear emission models 
that predict a strong dependence on the form of nuclear power out- 
put with nuclear accr etion rate (iFalcke. K ording & Markof3 l20()i 
Churazov et al .120051) . Here we examine the nuclear emission prop- 
erties in a sample of over 50 BCGs using the energy demands of 
X-ray cavities as a measure of mean nuclear accretion rate. 

We assume H = 70 km s _1 Mpc -1 , Q m = 0.3 and £l A = 0.7. 
All errors are \o unless otherwise noted. 



2 DATA REDUCTION AND ANALYSIS 
2.1 Sample selection 

We intend to explore the emergent properties of accreting black 
holes at the centres of clusters. Therefore we have selected objects 
with large X-ray cavities which we use to estimate the mean accre- 
tion rate of each object over the past 10 7 — 10 8 yr. Sources with a 
range of cavity powers, and thus accretion rates, were selected. The 
sources were selected from cluster, group and elliptical galaxy sam- 
ples which sho w evidence of AGN activity in the form of cavities in 



X-ray images feirzan et al. 2004; Raffertv et al. 2 0061 : 1 Allen et alJ 
2006: ICavagnolo et al. 2010; O' Sullivan et al. 2011a|). These ob- 



jects were supplemented with other recently discovered X-ray point 
sources in cavity systems (RXCJ0352.9+1941, RXCJ1459.4-1811, 
RXCJ1524.2-3154, RXCJ1558.3-1410, Zw348) and three noncav- 
ity systems (Zw2089, A2667, A611) each with a bright point 
source. 

We have also included three quasars taken from the 

literature for compariso n, H1821+643, IRAS 09 104+4 109 

and 3C186 dRussell et al. l20ld : Cavagnolo et al. submitted; 
ISiemiginowska et alJl2005L boiQh . These sources have very differ- 
ent spectral energy distributions (SEDs) from the low luminosity 
AGN dominating this sample. All three quasars have broad optical 
emission lines and bolometric luminosities of ~ 10 47 erg s _1 , 
far greater than the rest of the sample. Both HI 821+643 and 
3C186 have strong big blue bumps in the optical-UV band 
whereas IRAS 09 104+4 109 is a heavily obscurred quasar with 
most of its bolometric luminosity emerging in the IR. The big 
blue bump emission is usually interpreted as thermal emission 
from an accretion disk around the SMBH, which has then been 
re-radiated in the infrared in IRAS 09 104+4 109. This component 
is absent in radiatively inefficient low luminosity AGN (eg. 
IChiaberge. Caoetti & Macchettoll2005h . 

This sample is neither complete nor unbiased and we are 
mindful of this in our interpretati on and analysis. Our sam ple over- 
laps the 15 objects analysed by IMerloni & Heinzl d2007h and we 
have quadrupled the sample size, including a large number of upper 
limits. In total, 57 sources were selected covering a redshift range 
from the nearest Virgo ellipticals at a distance of only 17 Mpc to 
3C186 at z = 1-06 and a mass range from single elliptical galaxies 
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to rich clusters. We sample mean accretion rates from 2 x 10 -6 to 
0.6MrHdd for the first time. The full list of sources, excluding the 
three quasars, is shown in Table Q] 



2.2 Chandra data reduction 

For each object in this sample, we selected the deepest Chandra 
observation in the archive for analysis (Table Q]). Each observation 
was reprocessed using CIAO 4.4 and CALDB 4.4.7 provided by 
the Chandra X-ray Center (CXC). The level 1 event files were re- 
processed to apply the latest gain and charge transfer inefficiency 
correction and then filtered to remove photons detected with bad 
grades. The improved background screening provided by VFAINT 
mode was also applied where available. Background light curves 
were extracted from the level 2 event files of neighbouring chips 
for observations on ACIS-I and from ACIS-S1 for observations 
on ACIS-S3. The background light curves were filtered using the 
LC .CLEAN scripQ provided by M. Markevitch to identify periods 
affected by flares. The final cleaned exposure times of each ob- 
servation are detailed in Table Q] Standard blank- sky backgrounds 
were extracted for each observation, processed identically to the 
events file and reprojected to the corresponding sky position. The 
blank- sky background was then normalized to match the count rate 
in the 9.5 — 12 keV energy band in the observed dataset. This cor- 
rection was less than 10 per cent for the majority of the observa- 
tions. Each normalized blank sky background was also checked 
against the observed background spectrum extracted from a source- 
free region of each dataset to ensure it was a good match. 



2.3 X-ray point source flux 

Fig. Q] shows an example of a clear detection of surface brightness 
depressions ind icating cavities and an X-ray point source at the cen- 
tre of A2052 telanton. Sarazin & McNamarall2003l : iBlanton et all 
2009, |20T ]]). We required at least a 3cr detection above the 
background, where the background error considers only Pois- 
son statistics, in a hard 3 — 7 keV energy band image to con- 
firm the detection of an X-ray point source. The majority of 
the sources were not detected with sufficient counts above the 
cluster background to generate a reasonable spectrum. Therefore, 
we used two alternative methods, based on those described by 
iHlavacek-Larrondo & Fabianl J2Q1 lh , to calculate the X-ray flux of 
the confirmed sources and la upper limits on the non-detections. 

The photometric method sums the point source flux in a 
1 arcsec radius region using an exposure-corrected image, which 
was weighted by a spectral model for the point source. This re- 
gion was centred on the peak in the X-ray cluster emission if 
no point source was detected. The only exception was the Cen- 
taurus clus ter where we used the position of the radio source as 
the centre (iTavlor et al.ll2006h . For the majority of the sample, we 
used an absorbed powerlaw model PHABS(ZPHABS(POWERLAW)) 
with no intrinsic absorption, a photon index Y = 1.9 (eg. 



Gilli, Comastri & Hasinger 2007) and Galactic absorption from 



Kalberla et all (2005). For sources detected with several hundred 



the C-statistic available in XSPEC to determine the best-fit pa- 
rameters for spectra with a low number of counts dCashl [l979l ; 
Wacht er. Leach & Kell ogdl979h . Note that the intrinsic absorption 
was set to zero for sources with only upper limits and the pho- 
ton index was set to 1.9 where this was consistent with the best-fit 
value within the error. The cluster background was subtracted us- 
ing an annulus around the point source from 1.5 — 2.5 arcsec and an 
exposure-corrected image that was uniformly weighted at the aver- 
age peak in the cluster spectrum around 1.5 keV. The point source 
flux was also corrected for the fraction of the point spread function 
(PSF) falling in the 1 arcsec region ( 90 per cenQ). A larger region 
encompassing more of the PSF also included a greater fraction of 
cluster background, which increased the measurement uncertainty. 

However, for galaxy clusters and groups with steep central 
surface brightness peaks this method is likely to significantly un- 
dersubtract the background cluster emission. We therefore also em- 
ployed a spectroscopic method where we fit a spectrum extracted 
from the 1 arcsec radius point source region with a model for both 
the point source and cluster emission. The parameters for the clus- 
ter model were determined by extrapolating profiles of the pro- 
jected cluster properties in to the point source region. These pro- 
files were generated by extracting spectra from a series of circular 
annuli centred on the point source and fitting them with a suitable 
spectral model to determine the temperature, metallicity and nor- 
malization. We required a minimum of ~ 2000 source counts per 
region to ensure good constraints on the cluster parameters. How- 
ever, fewer counts per region were allowed for the low temperature 
sources in the sample where the Fe L line emission improves tem- 
perature diagnostics. Point sources were identified using the CIAO 
algorith m WAVDETECT, visua lly confirmed and excluded from the 
analysis dFreeman et al.ll2002h . Regions of nonthermal jet emission 
were also excluded. The cluster spectra were grouped to contain 
a minimum of 20 counts per spectral channel, restricted to the en- 
ergy range 0.5 — 7 keV and fit in XSPEC with appropriate responses, 
ancillary responses and backgrounds. 

For the majority of the c lusters, an absorbed single temper- 
ature PHABS (MEKAL) model (Balucinska-Church & McCammon 
1 992 ; iMewe. Gronenschild & van den Oordlll985l: ll986l~lKaastra 



1992; LiedahL Osterheld & Goldstein 1995) provided a good fit to 



source counts or more, this mod el was fitted t o an extracted spec- 
trum using XSPEC version 12 dArnaudl [19%! ) with the intrinsic 
absorption and the photon index left free. These model param- 
eters are detailed in Table Q] We used the modified version of 



the cluster emission in each annulus. This model was insufficient 
for some of the bright, nearby clusters with multiphase gas signa- 
tures, such as Centaurus, M87 and A2052. In these cases, a second 
MEKAL component or a MKCFLOW component was added to the 
spectral model. The cluster redshift and Galactic column density 
were fixed to the values given in Table [T] Abundances were mea - 
sured assuming the abundance ratios of lAnders & Grevessel dl989h . 
For the multi-temperature models, the metallicity was tied between 
the two components and, for the MKCFLOW model, the lower tem- 
perature was fixed to 0.1 keV and the higher temperature was tied 
to that of MEKAL component. We produced radial profiles for the 
best fit temperature, metallicity and normalization parameters and 
used powerlaw fits to the inner points to extrapolate these proper- 
ties into the point source region. The cluster spectra were also used 
to determine the Galactic absorption for sources located on regions 
of the sky where the absorption is highly variable. The spectral fits 
were repeated for annuli at large radii, excluding cooler gas com- 
ponents in the cluster core, leaving the nn parameter free (Table 

The spectrum extracted from the point source region was then 



1 See http://cxc.harvard.edu/contrib/maxim/acisbg/ 
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Figure 1. Exposure-corrected Chandra images covering the same field of Abell 2052 (see Blanton, Sarazin & McNamara 2003; Blanton et al. 2009, 2011). 
The colour bar has units photons cm -2 s -1 pixel -1 . Left: 0.5 - 7 keV energy band showing the X-ray cavities. Right: 3-7 keV energy band showing the 
AGN point source detection. 



fitted with a combined model for both the point source and cluster 
emission. The similarity between the powerlaw and thermal model 
components, particularly for higher temperature clusters, can make 
it difficult to distinguish them with Chandra's spectral resolution. 
Most of the model parameters were therefore constrained for the 
fit. The parameters for the point source model were set to those de- 
tailed in Table [T] with only the normalization left free. The cluster 
model temperature was fixed to the value determined by extrapo- 
lating from neighbouring annuli. The metallicity was fixed to the 
value in the neighbouring cluster annulus as it was generally not 
found to vary significantly at these small radii. The normalization 
of the cluster component was problematic because it was strongly 
affected by cavity substructure in the core, which produced large 
variation. We therefore constrained the cluster normalization in the 
point source region to be no less than the normalization from the 
neighbouring annulus, scaled by the ratio of their respective areas. 
The XSPEC CFLUX model was used to determine the flux or an up- 
per limit for the unabsorbed point source component (Table Q]). 

The main source of error in both the photometric and spec- 
troscopic measurements of the point source flux is the subtraction 
of the cluster emission. The photometric method is likely to under- 
estimate the background cluster emission and should therefore be 
treated as an upper limit on the point source flux. The spectroscopic 
method improves on this by allowing for the increase in cluster 
surface brightness towards the cluster centre but may significantly 
overestimate the cluster background because it is indistinguishable 
from the powerlaw component. For strongly obscured point sources 
and low temperature cluster emission, the spectroscopic method is 
likely to be a significant improvement over the photometric method. 
For higher temperature clusters, the photometric method may be 
more accurate. We have therefore listed both the photometric and 
spectroscopic fluxes in Table Q] but used the generally more accu- 
rate spectroscopic flux in our analysis. The possible bias in this 
measurement for higher temperature clusters with strong 3 — 7 keV 
emission is discussed in section [331 

Several of the brighter point sources in the sample were signif- 
icantly piled up in the longest Chandra exposures initially selected 
for analysis. Pile up occurs whenever two or more photons, arriving 
in the same detector region and within a single ACIS frame inte- 



gration time, are detected as a single event (Davis 2001). For M87 
and Cygnus A, there were alternative observations available in the 
archive with shorter 0.4 s frame times for which the point source 
was not piled up. These short frame time observations were used 
to calculate the point source flux and the cluster background was 
analysed using the deeper exposures. All the archival observations 
of the point source in Perseus, where the cluster centre is not posi- 
tioned far off axis distorting the PSF, were found to be significantly 
piled up. Perseus was therefore excluded from this sample. 

Fig. [21 compares our spectroscopic point source fluxes with 
measurements for the same sources available in the literature. The 
majority of the fluxes are consistent within the errors. Small varia- 
tions are expected due to differences in background subtraction and 
the position selected for the upper limits but there are three sources, 
Centaurus, NGC4782 and RBS797, with significantly different val- 
ues which we have considered in detail. There is only a modest dis- 
crepancy for RBS797 given the large errors and this is likely due 
to the additional model components used for the ICavagnolo et alj 
(l201ll) result. The differences for Centaurus and NGC4782 are 
due to our use of a two temperature rather than a single temper- 
ature cluster model. The best-fit single temperature falls midway 
between the preferred higher and lower temperature values of the 
two component model and therefore significantly underestimates 
the cluster surface brightness in the 2 — 10 keV band used to de- 
termine the point source flux. Our two component model therefore 
finds a higher cluster background and a significantly lower nuclear 
point source upper limit. 

2.4 Cavity power 

The cavities observed in the X-ray images of this sample allow 
a direct measurement of the mechanical ou t put from the AGN 
dChurazov et all l2000l: iDunn & Fabianl 120041 : lBirzanetal.1 12004 
McN amara & NulsenN2007h . For a bubble filled with relativistic 
plasma, the energy required to inflate it is given by E = 4PV, where 
the bubble is assumed to be in pressure equilibrium with the sur- 
rounding ICM. The bubble energy is then divided by the sound 
speed timescale or the bu oyant rise time to estimate the power in- 
put to the ICM (see eg. Birzan et al. 2004). Cavity powers for the 
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Figure 2. Comparison of the spectroscopic point source fluxes with values 
from the literature (see Table [!}• Note that the point source fluxes shown 
here were evaluated in the same energy band as those given in the literature 
and may differ significantly from those in Table [T] 



majority of the targets in this sample w ere available in the liter- 
ature and otherwise estimated following iBirzan etlril (2004) (see 
Table Q]). These values should be treated generally as lower lim- 
its on the total mechanical energy input by the central AGN. The 
datasets available in the Chandra archive for this cluster sample 
vary from snapshot 10 ks exposures to almost complete orbit ex- 
posures over 100 ks. Even for the nearest and brightest galaxy 
clusters, deeper ex posures continue to discover new cavities (eg. 
Fabian et al. 2011). Some cavities will therefore have been missed 
from these estimates of the total power. In addition, weak shocks 
and sound waves have also been found to contribute significantly to 
the power output of the central AGN ( eg. IFabian et al Jl2003l l2006: 
iForman et al"]l2005h . 

Estimates of the cavity power can be complicated by the pres- 
ence of significant X-ray emission from jets. For example, the syn- 
chrotron self-Compton emission in the centre of 3C 295 makes it 
difficult to determine the extent of the inner cavities and there- 
fore the total power is significantly underestimated for this sys- 
tem (lHarris et alJl200d : lAllen. Ettori & Fabianll200ll) . The cavity 
power estimate for 3C295 is therefore shown as a lower limit. 
Three sources in our sample, Zw2089, A611 and A2667, have 
point sources but no detected cavities in the X-ray images. Ideally, 
these sources would be included in the analysis with upper limits 
on the cavity power. However, it is difficult to place meaningful 
constraints on the possible size of non-detected cavities, particu- 
larly at large radius where the noise increases and the X-ray surface 
brightness drops off rapidly (see Birzan et al. 2012). These sources 
have therefore been included with only illustrative upper limits, not 
quantitative. 

By comparing the radiative output from the X-ray point source 
with the cavity power, we can calculate the radiative efficiency of 
the AGN. The total power output from the AGN can also be used 
to infer a mean accretion rate, given an assumption of the accre- 
tion efficiency. We refer throughout to a mean rather than instanta- 



neous accretion rate as this method estimates the average accretion 
requirements over the ~ 10 7 yr age of each cavity. We have also 
calculated the theoretical Eddington and Bondi accretion rates for 
comparison with the inferred accretion rates for the sample. 



2.5 Eddington luminosity 

The Eddington luminosity indicates the limiting luminosity of the 
SMBH when the outward pressure of radiation prevents the gravi- 
tational infall of accreting material. For a fully ionized plasma, the 
Eddington luminosity can be expressed as 



=1.26xl0 47 



ergs 



10 9 M G 



(D 



where Mbh is the SMBH mass. Dynamical estimates of 
the black h ole mass were used where these were av ailable 
(Cygnus A, [f adhunter etaO 120031: iRaffertv et all 120061 : M84, 



Walsh, Barth & Sarzi 2010; M87, Gebhardt et al. 2011; NGC4261 



Ferrarese, Ford & Jaffe 1996). For the majority of the sample with- 



out dynamical masses, we relied on well-studied relations between 
the b lack ho l e mass and the properties of the host galaxy. Follow- 
ing iGrahaml d2007h . the apparent K-ba nd magnitude of the host 
galaxy from the 2MASS catalogue^ dSkrutskie et all [2006) was 
converted to an absolute magnitude, correc ted for Galactic ex- 
tinction (Schlegel, Finkb emer & Davis1 [l998), redshift and evolu- 
tion ( Poggiantilll9970 . and used to estimate the black hole mass. 
3C 401, RBS797 and Zw 348 have detected X-ray point sources but 
no available apparent K-band magnitudes for this analysis and are 
not incl uded in plots with Edd ington accretion rate or luminosity 
scaling. iDalla Bonta et al.l (120091) find an upper limit on the central 
black hole mass in A2052 from dynamical measurements, which 
is consistent with the value calculated from the apparent K-band 
magnitude. 

There are however potential problems with the use of K- 
band magnitudes to estimate black hole masses. ILauer et al 
suggest that the apparent magnitudes from 2MASS are not deep 
enough to capture the extent of the BCG envelope and there- 
fore will underestimate the total luminosity and black hole mass. 
iBatcheldor et al.l d2007h instead suggest that the BCG extended en- 
velope, formed of debris from tidal stripping, is unlikely to be 
closely associated with the central galaxy dynamics. We have there- 
fore considered these black hole masses to be estimates and, where 
possible, present the results with and without scaling by the Ed- 
dington accretion rate. 



2.6 Bondi accretion rate 

The primary aim of this paper is to examine the relationship be- 
tween accretion and power output in SMBHs. One of the outstand- 
ing problems is whether the jets are powered by cold accretion 
from circumnuclear accretion disks of atomic and molecular gas 
or whether they are fuelled by a nearly spherical inflow of hot keV 
gas. It is impossible to prove either case because we lack imaging 
on the scale of and below the Bondi radius for almost all objects in 
our sample. However this has been done in nearby ellipticals where 
Chandra's resolution is close to the size of the Bondi sphere (eg. 
| Di Matteo etaflbOOlL 120031 : IChurazov et ai1l2002l: IPellegrini et al.1 
2003). This problem has been addressed by I Allen et al.l (2006) in a 
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sample of nearby ellipticals and bv lRaffertv etail d2006h m a more 
distant sample with properties similar to our o wn and these stud - 
ies arrived at somewhat different conclusions. lAUen et alJ d2006h 
found a strong trend between jet power and an estimate of the Bondi 
accretion rate based on circumnuclear X- ray properties whereas the 
sample analysed by Raffe rtv et al.l {2006) required greater extrapo- 
lation in to the Bondi sphere making their conclusions on Bondi ac- 
cretion highly uncertain. Nevertheless, using energetic arguments, 
they found Bondi power accretion was unable to fuel the most pow- 
erful jets unless their black holes were much larger than implied by 
the M-cj relation. This question is crucial for understanding the re- 
lationship between jet power and nuclear emission therefore here 
we re-analyse data for a subsample of 13 systems where we have 
resolution close to the Bondi sphere. Thi s subsample includes the 
9 galaxies targeted by All en et alJ (2006) and 4 additional objects 
with new, deeper Chandra observations (Table E}. 

Assuming spherical symmetry and negligible angular momen- 
tum, the Bondi rate, Mb, is the accretion rate for a black hol e em- 
bedde d in an atmosphere of temperature, T, and density, n e , teondil 
Il9521) and can be expressed as 

M yr- 1 V keV; V cm" 3 J \ 1O 9 M J 

for an adiabatic index 7 = 5/3. This accretion occurs within the 
Bondi radius, tb, where the gravitational potential of the black hole 
dominates over the thermal energy of the surrounding gas, 



temperature profiles were generally found to flatten in the central 
regions of these systems. We have therefore assumed that the tem- 
perature at the accretion radius has only decreased by a further 
factor of 2 from the innermost temperature bin. Using equation [3] 
the Bondi radius was calculated from this innermost temperature 
value and the black hole mass. Dynamical black hole mass esti- 
mates were used where available. For comparison with I Allen et al. 
( 2006), stellar velocity dispersions from the HyperLeda Databas^ 
were then used to calculate black hole masses for the rest of the 
subsample (iTremaine et al]|2002h . Only HCG62 did not have a dy- 
namical or velocity mass estimate and th e K band magnitude was 
used (section 12.51 ). ITremaine et al] d2002h find an intrinsic disper- 
sion of 0.25 — 0.3 dex in log(MBH) in the Mbh — o relation, which 
dominates over the error in the measurement of the velocity disper- 
sion. 

The deprojected electron density continues to increase in the 
cluster centre and was therefore extrapolated to the Bondi radius 
using several different model profiles. Three different models were 
considered: a powerlaw model continuing a steep density gradient 
to t*b, a P -model flattening to a constant and a shallowing Sersic 
profile with n = 4. These models were fitted to the density profile 
and used to calculate the gas density at r%. The density at r% is 
shown as a range of likely values from the powerlaw model upper 
limit to the /3 -model or Sersic model lower limit. Using equation 
13 the Bondi accretion rate was calculated from the density and 
temperature at r% and the black hole mass. 



kpc \kcVj \W>M e J- W 

The Bondi accretion rate is therefore an estimate of the rate of 
accretion from the hot ICM directly onto the black hole and de- 
pends on the temperature and density of the cluster atmosphere at 
the Bondi radius. 

The projected cluster spectra, extracted from a series of 
annuli as described in section 12.31 were deprojected using 
the model-independent spectral deprojection routine DSD EPROJ 
dSanders & Fabianl l200l iRussell. Sanders & Fabianl l2008h . As- 
suming only spherical symmetry, DSDEPROJ starts from the 
background- subtracted spectra and uses a geo metric method 
(iFabian et al.lll98ll : lKriss. Cioffi & Canizaresll 19831) to subtract the 
projected emission off the spectrum from each successive annulus. 
The resulting deprojected spectra were each fitted in XSPEC with an 
absorbed single temperature MEKAL model to determine the tem- 
perature and density of the gas, as described in section [231 Several 
of the selected clusters have clear evidence f or multiple tempera- 
ture comp onents in the inner r egions (eg. M87 jForman et al.ll 2005; 
Centaurus jFabian et al.ll2005h and it is likely that deeper exposures 
of other clusters in the subsample will also produce robust detec- 
tions of multi-phase gas. However, given the available range in ex- 
posure depth, a uniform method with a single temperature model 
was used to determine the emission-weighted average temperature 
of the ICM in the cluster centre. We also generated deprojected 
electron density profiles with finer radial binning from the surface 
brightness profile s and incorporating the temperature and metallic- 
ity variations (eg. lCavagnolo et aDl2Q09h . 

These profiles trace the cluster parameters to radii within an 
order of magnitude of the Bondi radius and therefore some extrap- 
olation is required. Note that although the Bondi radius in M87 
is resolved by Chandra, sl significant region is affected by pileup 
and the PSF from the jet knot HST-1 and must be excluded. The 



2.7 Radio point source flux 

We also compared the nuclear X-ray flux with the radio core flux 
for a subset of the sample to try to determine the origin of the 
nuclear X-ray emission. Radio observations were available for 22 
sources in our sample which allowed us to reliably distinguish on- 
going core activity. If the sole source of the X-ray flux is in the base 
of a jet, a direct relationship between the radio and X-ray core flux 
would be expected. 

Nine of the sources in this subsample were observed simulta- 
neously at C and X bands with the ATCA (project C1958, PI Edge). 
All but one of the remaining sources were observed at C-band with 
the VLA-C array (various projects, PI Edge), with the last source 
(A478) having been observed simultaneously at L and X band with 
the VLA-A array (project AE117). For each BCG, the radio-SED 
was further populated with data from the major radio catalogs (in- 
cluding but not limited to: AT20G at 20 GHz, NVSS/FIRST at 
1.4 GHz, SUMSS at 0.843 GHz, WENSS/WISH at 0.325 GHz, 
VLSS at 0.074 GHz). Additional fluxes were found by searches 
around the radio-peak coordinates in both the NED and HEAS ARC 
online databases. All literature fluxes were individually scrutinised 
to ensure reliable matches. Where the synthesised beam size was 
considered limiting, leading to source confusion, data were dis- 
carded. Four of the sources have VLBA observations at C-band. 

Core flux contributions were calculated by considering both 
the morphology and SEDs of each of the sources. The VLBA C- 
band observations provided direct measurements of the core flux. 
For the remainder, the SEDs were decomposed into two major com- 
ponents; a flatter spectrum, active component attributed to ongo- 
ing activity within the core of the AGN and a steeper spectrum 
component, most dominant at lower frequencies attributed to either 
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Figure 3. Histogram showing the distribution of the ratio of X-ray point 
source luminosity to cavity power for the low luminosity AGN in the sam- 
ple. Three point sources and one upper limit have no detected cavities and 
are therefore not included. 



past AGN activity or alternate acceleration mechanisms (e.g. radio 
lobes, mini-haloes, etc.). 

Where a clearly resolved core was present in the observations, 
two component SEDs were fitted directly. For sources which were 
resolution limited at at C-band, spectral breakdown of the SEDs 
was performed on a case-by-case basis. Consideration was given 
to extent seen at other wavelengths, spectral shape and variabil- 
ity, with the proviso that variable sources are more likely contain 
a strong currently active core. Simple mathematical models were 
fit to the SEDs using IDL routines where a strong case could be 
made for believing distinct components were present. For sources 
where past and current activity could not be reliably distinguished, 
limits were placed on the core contribution. Full details of the SED 
analysis will be presented in Hogan et al. (in prep). 



3 RESULTS 

In total, 27 out of 54 BCGs in this sample were found to have X- 
ray central point sources detected above 3 a in the 3 — 7 keV energy 
band. Although this is not a complete sample of objects, central 
X-ray point sources appear to be common in BCGs with detected 
X-ray cavities. Fig. [3] shows the distribution of the ratio of X-ray 
nuclear luminosity in the 2 — 10 keV energy range to cavity power, 
Lx/Pcav, for the cluster sample. The detected point sources in this 
sample cover a broad range in radiative efficiency from Lx/P C av > 
0.1 (eg. 3C295) to L x /P C av < 10~ 3 (eg. A2199). This distribution 
presumably reflects a broad range in average accretion rate with 
the majority accreting in a radiatively inefficient mode ie. they are 
ADAFs (see section [T]). 



3.1 Radiative and cavity power output 

Fig. |4] (left) shows a correlation between the nuclear point 
source 2 — 10 keV luminosity and the cavity power injected 



into the surround ing ICM. This correlation was first found by 
Merl oni & Heinzl d2007h for a sample of 15 AGN with measured 
cavity powers, 13 of which also had X-ray point source detections. 
The scatter in this correlation covers three orders of magnitude 
therefore the gen eralised Kendall's t rank correla tion coefficient 
for censored data dBrown. Hollander & Korwarll974 ) from the sur- 
vival analysis package ASURV Rev 1.3 dlsobe. Feig elson & Nelson 
1986; Isobe & Feigelson 1990; Lavalley, Isobe & Feigelson 1992) 
was used to evaluate its significance. We find a probability of ac- 
cepting the null hypothesis that there is no correlation between the 
point source luminosity and the cavity power of P nu \\ = 4 x 10 -4 . 
Given the substantial difference in the timescales, ie. six orders of 
magnitude, it is surprising to observe a trend at all between these 
properties and suggests that AGN feedback is persistant. 

It is not clear if the systems with upper limits on the point 
source luminosity form part of this trend, and are only just too 
faint to detect, or are currently 'off. We have therefore consid- 
ered these two scenarios separately using the BCES estimators 
dAkritas & B ershadv 19961) for a linear regression fit to the data 
points only and the ASURV non-parametric Bu ckley- James linear 
regre ssion fit to both the points and upper limits dBucklev & Jamesl 
119791) . Note that this analysis excluded the point sources with no 
detected cavities in the surrounding ICM as no effective upper lim- 
its for the cavity power could be estimated (see section 12.41 ). The 
inclusion of the upper limits produces a significant shift in the cor- 
relation but only a small change in the observed slope, which is 
consistent within the error. The BCES orthogonal best fit was found 
to be log(P C av) = (0.69 ± 0.08)log(L x ) + 15 ± 3 compared to the 
ASURV Buckley- James best fit of log(P cav ) = (0.7 ± 0. l)log(L x ) + 
14. The majority of the sample analysed are in a radiatively ineffi- 
cient accretion mode where the mechanical cavity power dominates 
over the radiative output from the AGN. The observed slope shows 
that the radiative efficiency of the X-ray nucleus increases with in- 
creasing cavity power. 

For comparison, we have included three quasars from the lit- 
erature that are located in luminous cool core clu sters to illus- 
trate t his increase in radiative efficiency. HI 82 1+643 dRussell et al.l 
l2010h. IRAS 09104+4109 (Cavag nolo et al. submitted) and 3C 186 
dSiemiginowska et al.l2005L[20101) each have a total radiative power 
of L ~ 10 47 erg s _1 that exceeds their cavity power by at least an 
order of magnitude. Fig.|4](left) shows that these sources appear to 
form an extension of the trend to much higher cavity powers and 
presumably much higher accretion rates. This implies that as the 
accretion rate r ises black holes become more radiatively efficient. 

Following Mer loni & Heinzl d2007h . we also estimate the bolo- 
metric point source luminosity for all sources in the sample and 
scale the luminosity and cavity power by the Eddington luminosity 
(Fig.|4]right). The SED for low luminosity AGN lacks the 'big blue 
bump' of emission dominating higher accretion rate sources and 
is likely to be dominated b y the emission at hard X-ray energies. 
IVasudevan & Fa bian (2007) find a typical bolometric correction for 
low luminosity sources of order ~ 10, however for e ase of compar- 
ison we adopt the same value as Merloni & Heinz] d2007h . Fig. ED 
(right) shows the correlation between the bolometric point source 
luminosity and the cavity power where both quantities are scaled 
by the Eddington luminosity (section 12.51 ). The non-parametric 
Buckley-James linear regression method from ASURV software 
was used to deter mine the best-fit relatio n for both the points and 
the upper limits dBucklev & Jamesl 1 19791) . The best-fit slope of 
log(L cav /L Edd ) oc 0.91og(L nuc /L R dd) is significantly steeper than 
that found by iMerloni & Heinzl d2007h (Fig. @] right). Using the 
BCES estimators linear regression fit to only the data points we 
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Nuclear luminosity 2- 1 OkeV (erg/s) log(LMjc/L EDD ) 

Figure 4. Left: Nuclear X-ray luminosity in the energy range 2 — 10 keV calculated using the spectroscopic method versus cavity power from the literature. 
Right: bolometric nuclear luminosity versus cavity power where both quantities are scaled by the Eddington luminosity. Clusters with confirmed point source 
detections are shown by the filled circles and upper limits are shown by the open circles. Three quasar sources from the literature are includ ed for comparison 
(red cr osses). The variable sources are shown as blue stars (see section [331 The best-fit to both points and upper limits using the ASURV Buc kley & Jamesl 
( 1979) estimator for cen sored data is shown by the solid line. The BCES orthogonal fit to the data points only is shown by the dashed line. The best-fit from 
Merloni & Heinz ( 2007) is shown as a dotted line and is consistent with the BCES fit to the data points within the errors. 



determ ine that the best fit is consistent with the iMerloni & Heinzl 
d2007h result within the errors. The inclusion of a large number 
of upper limits results in a significant difference in the slope of 
the correlation. We show in section [331 that 3 systems vary by up 
to an order of magnitude over a ten year timespan. This is appar- 
ently a significan t ly con tributing factor to the scatter to which the 
Merl oni & Heinzl J2007h sample was not sensitive. Another factor 
that may be contributing is b eaming of the central X -ray source, 
which is discussed in detail in Merloni & Heinz (2003). 

3.2 Selection effects 

It is also clear that the scatter in this correlation between the nuclear 
X-ray luminosity and the cavity power is underestimated by our 
sample selection. Sources were selected primarily on the detection 
of cavities in X-ray observations and therefore sources with point 
sources but no cavities are generally missed from the lower right of 
Fig- HI Three such sources, Zw 2089, A61 1 and A2667, were added 
to the sample to illustrate this selection bias. Note that the upper 
limits on the cavity powers are only illustrative, not quantitative 
(section |24l ). Whilst the observation of A2667 is shallow at only 
8 ks, both A611 and Zw2089 have sufficiently deep X-ray obser- 
vations that would detect signs of any cavities or feedback-related 
substructure in the cluster cores. In particular, Zw 2089 contains a 
bright X-ray point source and appears very relaxed with smooth 
extended X-ray emission and no likely cavity structures. However, 
if cavities in this system are emerging along our line of sight they 
will be particularly difficult to detect. There could be a significant 
number of similar systems with bright point sources but no cavity 
structures, which will tend to increase the scatter in the observed 
correlation further. 

The brightness and temperature of the surrounding cluster 



emission could also introduce another selection effect to this anal- 
ysis. Point sources are identified by a significant detection of emis- 
sion in a hard X-ray band, 3 — 7 keV, above the background clus- 
ter (section |231 >. Bright, high temperature clusters will have more 
emission in this energy band than fainter, cooler systems therefore 
potentially making it more difficult to detect a point source above 
this background. Fig. [5] shows that for systems with strong back- 
ground cluster emission in the 3 — 7 keV band there are several 
sources (labelled) with higher upper limits. This suggests that an 
AGN in these BCGs would have to be brighter to be detected than 
in BCGs with fainter emission in this energy band. However, only 
a handful of sources in our sample appear to be affected so this will 
only slightly reduce the scatter in Fig.|4] 

The selection bias in our sample will have a significant im- 
pact on the best-fit linear relation determined for the correlation in 
Fig. [4] (right). This is illustrated by the difference in slope between 
the Buckl ey- James linear regressio n in this analysis and the result 
found by Merl oni & Heinzl d2007h . We have therefore not drawn 
any further conclusions from the slope of this correlation but con- 
sidered the possible sources of the scatter, which covers at least 
three orders of magnitude. 

3.3 Uncertainty in cavity power and black hole masses 

The measured cavity powers available in the literature for a par- 
ticular system are often found to differ by factors of a few up to 
an order of magnitude. This can be due to new observations of an 
object, which reveal more cavities or better constrain the shapes of 
previously known cavities. It also reflects the inherent systematic 
un certainty and judgement of t he extent of the cavity volume (see 
eg. iMcNamara & Nul sen 2012 ). Cavities with bri ght rims, such as 
those in the Hydra A (eg. IMcNamara et alJ[2QQ0h and A2052 (eg. 
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Figure 5. Nuclear X-ray flux in the 2—10 keV energy band versus the 
cluster flux from a surrounding annulus in the 3 — 7 keV energy band. 



Blantonetal. 2011), have a well-defined shape, although some un- 
certainty still exists over the line of sight extent and whether the in- 
ner, outer or middle of the rims should be used. Most cavities do not 
have complete rims and their extent is difficult to constrain given 
the rapid decline in X-ray surface brightness with radius (Birzan et 
al. 2012). The cavity power for A2390 has a particularly large error, 



Pc?n = 1-0-to 9 x erg s _i , because the extent of cavities is dif- 
ficult to determine from the X-ray image. As discussed in section 
12.41 the cavity power for 3C295 is likely to have been significantly 
underestimated and this is therefore an outlier. Therefore, for the 
majority of the sources, uncertainty in the cavity power is unlikely 
to introduce scatter greater than an order of magnitude. Also, al- 
though the black hole masses are likely to be a significant source 
of additional error for Fig. |4] (right), the scatter in Fig. |4] (left) is 
comparably large and does not depend on black hole mass. 



3.4 Absorption 

Another significant source of uncertainty in the X-ray point source 
fluxes may be attributable to photoelectric absorption by interven- 
ing gas within the galaxy or in a circumnuclear torus. This would 
cause a systematic underestimate of the point source luminosity. 
The amount of intrinsic absorption can be determined from spectral 
fitting but only for the brighter sources in our sample. Fig. [6l shows 
the intrinsic absorbing column density as a function of nuclear X- 
ray luminosity for 25 of the 27 detected low luminosity AGN. There 
were insufficient counts for the detections of 3C388 and A2667 to 
constrain the intrinsic absorption. The intrinsic nn values cover a 
range from relatively unobscured sources with nn < 10 21 cm -2 , 
such as M87, up to heavily obscured narrow line radio galaxies 
with nu > 10 23 cm -2 , such as 3C 295. It is therefore plausible that 
a number of the non-detected sources could be moderately or heav- 
ily absorbed. This would account for some of the scatter in Fig. |4] 
and potentially render these sources undetectable in the 0.5 — 7 keV 
energy band accessible to Chandra. 

Fig. [6] has interesting astrophysical implications as well. The 



Figure 6. Nuclear X-ray luminosity versus the intrinsic absorption of the 
point source for a subset of objects with sufficient counts for a reasonable 
spectral fit. 



most luminous sources in our sample often contain large columns 
of cold intervening gas along the line of sight, in some cases 
exceeding 1 23 atoms cm -2 . It is unknown where along the line 
of sight this gas lies but it is likely to be close to the nu- 
cleus. This implies the existence of a supply of cold gas to 
fuel the nucleus. The range of intrinsic absorption may be re- 
lated in part to the geometry of cold gas relative to the cen- 
tral point source. For example the point sources in Hydra A and 
NGC4261 are peering through circumnuclea r gas disks that are 
highly inclined to the plane of the sky (eg. | jaffe & McNamaral 



ll994lDwarakanath. Owen & van Gorkoml 1995k Furthermore the 
path towards the core in Cygnus A is strongly reddened by inter- 
vening dust and presu mably accompanying gas loca ted within 800 
par sees of the nucleus dVestergaard & Barthelll 19931) . 

Using a subset of example upper limits at different redshifts, 
we determined that if these undetected sources had an intrinsic ab- 
sorption of 10 22 — 10 23 cm -2 their true luminosity could be fac- 
tors of up to a few greater than the observed upper limit. Intrinsic 
absorption of ~ 10 24 cm -2 and above is required for the true lu- 
minosity of the source to be an order of magnitude or more greater 
than the observed upper limit. For a significant amount of the scat- 
ter in Fig.|4]to be generated by intrinsic absorption, a large fraction 
of the upper limits must therefore be Compton thick. These sources 
would then represent a very different population to the detected 
point sources in this sample. 

Although a large fraction of the upper limits in this 
samp l e could be significantly absorbed (eg. Maiolino _et al.1 
I l998l: [Risa liti. Maiol ino & Salvatilll999l: iFabian & Iwasawalll999l : 
Brandt & Hasi ngeJ2005l : lGuainazzL Matt & Perolal2005l) , this still 
seems unlikely to completely account for the three to four orders 
of magnitude scatter in Fig. |4] Reliably identifying Compton-thick 
AGN and dete r mining their intrinsic luminosity is difficult (eg . 
IComa stri 2004: iNandra & Iwasawall2007l: [Alexander et all [2 008). 
iGandhi et al.1 (12009) (see also lHardcastle. Evans & Crostonl 12009) 
found a linear relation between the 2—10 keV X-ray luminosity of 
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the AGN and the 12 jum mid -infrared luminosity. Somewhat sur- 
prisingly, [GancMeU^ J2009h also found that the 8 Compton-thick 
sources in their sample did not deviate significantly from this trend. 
We might therefore expect that a large fraction of the undetected 
sources in our sample have high infrared luminosities as the ab- 
sorbed power from the nucleus is re-radiated at these energies. 

However, for BCGs at the centre of cool core clusters, a 
significant fraction of the observed infrared luminosity is likely 
due to star format i on (eg. lEgami et alJl200d : lOuillen et aD 12008: 
lO'Dea etal.ll2008h . IE garni etalJ (120061) found that the BCGs in 
A1835, A2390 and Zw3146, all with point source upper limits 
in our sample, ar e infrared bright and have SEDs typical of star- 
forming galaxies. lOuillen et alj J2008h identify several BCGs as 
having very strong AGN contributions in the infrared, such as 
Zw2089, but in more modest cases this is difficult to disentan- 
gle from star formation. Interesting! v JOuillen et al.l d2008h find that 
the BCG in Sersic 159-03 is infrared faint with only an upper limit 
on the IR star formation rate. Therefore, despite being one of the 
strongest upper limits in our sample (Fig.©, Sersic 159-03 does not 
appear to host a heavily absorbed AGN. Without detailed SEDs to 
disentangle the AGN contribution from star formation in these ob- 
jects it is difficult to systematically determine if a significant frac- 
tion of the point source upper limits are heavily absorbed sources. 
However, even if this is the case, it is unlikel y to explain the th ree 
orders of magnitude scatter in Fig. HI (see also lEvans et al.ll2006h . 

3.5 X-ray variability 

Another source of scatter is variability in the nuclear X-ray lumi- 
nosity. The cavity power is averaged over the cavity ages, which 
are estimated from the sound cros sing time or the b uoyant rise 
time and are typically 10 7 — 10 8 yr (iBirzan et al]|2004l) . However, 
the nuclear power is expected to be variable on timescales much 
shorter than this. The effect of variability on th e trends in Fig. |4] 
is very much like that of rela tivistic beaming dMerloni & Heina 
l2007h . iMerloni & Heinzl (120071) further pointed out that variabil- 
ity is likely only to be a problem in the most luminous AGN such 
as Seyferts and quasars, which can vary on timescales of weeks. 
A subset of the clusters in our sample have multiple observations 
spaced by several years in the Chandra archive, which we searched 
for possible variations in the X-ray point source flux. 

The point source flux was calculated for each observation us- 
ing both the photometric and spectroscopic methods described in 
section |231 For the spectroscopic method, the differing depth of the 
observations could produce variations in the cluster parameters de- 
termined from surrounding annuli, which would add scatter to the 
measured point source fluxes. However, by fixing the cluster pa- 
rameters in each observation to those determined from the deepest 
exposure available, we verified that this did not significantly alter 
the results. The point source fluxes for this subsample are shown 
in Table [3] The best-fit intrinsic absorption and photon index are 
shown for each exposure but the fluxes were calculated using the 
values of these parameters from the deepest exposure of that source. 
No significant change in the intrinsic absorption or the photon in- 
dex was found for the sources analysed. 

Fig. [7] shows that the point source flux was found to signifi- 
cantly vary in A2052, Hydra A and M84. The flux in A2052 was 
observed to decline by an order of magnitude over the ten years 
traced by the Chandra archive. Hydra A shows a more modest de- 
cline by a factor of ~ 2 over 5 years, whereas for M84 a decline 
by a factor of ~ 3.5 drop in flux occurs in only 6 months. Note 
that for Hydra A the earliest ACIS-I observation in 2000 was not 



included (obs. ID 575) as this dataset was taken during the soft 
proton damage to the detector. This analysis is generally limited by 
the availability of suitably spaced observations of sufficient depth 
in the Chandra archive but suggests that a significant fraction of 
AGN in BCGs may be varying on timescales of months to a few 
years. The bright central point sources in the Perseus cluster and 
M87 have long been known to be variable at X-ray w avelengths 
(eg.|Rothschildetal. 1981: lHarris. Biretta & Ju nor 1997). Another 
source in our sample, NGC4261, has been found to be varia ble on 
short 3 — 5 ks timescales in a study bv lSambruna et alj (120031) . This 
variability is likely to be a significant source of scatter in the cor- 
relation between nuclear X-ray luminosity and cavity power. The 
cavity power is an average of the AGN activity over 10 7 — 10 8 yr 
whereas the point source luminosity is likely to fluctuate signifi- 
cantly on much shorter timescales potentially by orders of magni- 
tude. 

The shape of the nuclear spectrum from Hydra A is dramat- 
ically different from those of A2052 and M84. Hydra As spec- 
trum falls sharply below 2 keV while the flux below 2 keV in 
A2052 and M84 continues to rise. This strong decline in flux short- 
ward of 2 keV in Hydra A is due to a large column of intervening 
gas that may be associated with the large circumnuclear disk (eg. 
iDwarakanattu Owen & van Gorkomlll995l ; Hamer et al. in prep). 
We also searched for a change in the shape of the nuclear spectrum 
as the sources varied. Fig. [8] shows the spectrum for each observa- 
tion of each point source found to have significantly varying X-ray 
flux. The cluster background was subtracted from each spectrum 
using a surrounding annulus. The spectra are remarkably consistent 
between the observations and suggest that despite the large varia- 
tions in flux, particularly in A2052, there has been no significant 
change in the shape of the spectrum. 



3.6 Nuclear radio luminosity 

Fig. [9] (left) shows no apparent correlation between the nuclear X- 
ray flux and the 5 GHz radio core flux. There does appear to be an 
approximately linear trend between the nuclear X-ray luminosity 
and the radio luminosity (Fig. O right), although the X-ray flux is 
on average an order of magnitude larger. However, it is highly likely 
that this trend is due to redshift selection effects given the lack of a 
correlation in the flux-flux plot. 

Whilst care was taken to provide reliable core contributions to 
the overall radio flux density at C band, there are of course limita- 
tions. A variety of facilities were used to obtain the flux measure- 
ments used in the SEDs. Whilst this variety was considered in the 
decompositions, there will undoubtedly be situations where the true 
core contribution is lower than found in this analysis. This is due 
to contamination from extended emission in the lower resolution 
observations which is not adequately accounted for in the models. 
Similarly, for the highly core-dominated sources, large observed 
variability may lead to the radio core flux being underestimated at 
the epoch of the X-ray observations. These shortcomings will be a 
contributing factor to the scatter seen in Fig. [9] It should be noted 
however that the radio core contributions used here are taken from 
a larger sample of radio-loud BCGs analysed by Hogan et al (in 
prep). Of this larger sample, 26 are observed with the VLB A and 
strong agreement is seen between the direct VLBA core measure- 
ments and the SED-breakdown derived core contributions. Finally, 
many of the radio cores are self-absorbed so the 5 GHz flux may 
significantly underestimate the total radio power of the core. There 
is also likely to be significant scatter due to variability in both the 
X-ray and the radio flux. With no clear trend between the X-ray 
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Figure 7. AGN flux variability in A2052 (left), Hydra A (centre) and M84 (right). The point source fluxes are calculated using both the photometric (open 
crosses) and spectroscopic methods (filled circles). Note that a spectroscopic flux measurement could not be produced for M84 obs. ID 401 because the 1 ks 
exposure was too short. 




and radio nuclear flux it appears less likely that the X-ray emission 
originates solely from the base of a jet. 



3.7 Bondi accretion 

The deprojected temperature and electron density profiles for the 
Bondi subsample of 13 systems are shown in Fig. [TO] For each 
cluster we have marked the location of the Bondi radius and shown 
that the radial profiles are within roughly an order of magnitude of 
this. The two methods of calculating the deprojected density pro- 
file ^e_£onsisten^_as expected. For the clusters that overlap with 
the lAllen et all d2QQ6h sample, we generally find good agreement 
between the density and temperature profiles. The Cent aurus clus- 
ter pr ofiles were found to differ significantly because lAllen et all 
(2006) used a 35° wide sector to the NE of the nucleus com- 
pared to our full annuli, which included the complex structure W 
of the nucleus. The temperature profile for NGC5846 is also signif- 
icantly different in shape but we note that the central values used 
for the Bondi analysis are consistent. We used a more recent, deep 
observation of this sou rce and the results are consistent with the 
iMachacek et alj (1201 ll) analysis. 

Inner cavity substructure produced some sharp decreases 



in the deprojected density profile in several clusters, including 
NGC4636 and NGC5044 (Fig.[l0). The density models were there- 
fore fitted to all points within the central few kpc to smooth over 
substructure that is difficult to correctly deproject. In general, the 
inner radii of the density profiles were well-described by the three 
models used to extrapolate to the Bondi radius. The Bondi radius, 
accretion rate and cavity powers calculated from the temperature 
and density profiles f or each of the sel ected systems are shown 
in Table |2] Following lAllen et aD ^2006), we calculated the cav- 
ity power for only the inner two cavities of each object that are 
currently being inflated by the central AGN. 

Fig. Qj] (left) shows the Bondi power plotted against the cav- 
ity power generated by the inner two cavities in each s ystem. This 
show s a significant weakening of the trend found by lAllen et all 
( 2006) driven mainly by a difference in the estimates of cavity vol- 
ume. For most of these sources, our cavity powers are consistent 
with lAllen et alj d2006h within the large errors of a factor of 2 — 3 
on these values. However, for M84, M89, NGC4472 an d NGC507 
our cavity powers are lower than those of I Allen et al. (120061) by 



Merloni & Hein: 



factor s of up to an order of magnitude (see also 
l2007h . This was partly due to the availability of new, deeper ob- 
servations of M84, NGC507 and NGC4472, which more clearly 
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Figure 9. X-ray 2—10 keV vs. radio 5 GHz flux (left) or luminosity (right). Nuclear X-ray detections are shown by the solid points and X-ray upper limits are 
shown by the open points. 



revealed the cavity extent. lAUen etal.l(l2006h also used 1 .4 GHz ra- 
dio images to determine the edges of the cavities and this may have 
caused significant differences from our primarily X-ray method. 
Our new estimates of X - ray cavity pow er agree with esti mates from 
ICavagnolo eUulfcOlOh . lRaffertv et alJ d2QQ6h and b' Sullivan et al] 
d2011ah . It is not clear whether cavity powers will be more ac- 
curate when calculated using the X-ray or the radio observations, 
therefore we have also included Fig. Qj] (right) showing our anal- 
ysis of the Bondi accretion power versus the cavity po wers from 
Allen eTal1 ( [2006h . This plot shows a larger scatter than Alle n et aD 
(2006) found an d this scatter is due solely to differences in how we 
and I Allen et alJ d2006h calculated the central density. The vertical 
'error bars' should not be interpreted as such. Instead they repre- 
sent the range of Bondi powers from the three best-fit models and 
the midpoint is marked as no model provides a significantly bet- 
ter fit. The exception is the Centaurus cluster where the /3 -model 
is significantly preferred over the powerlaw and Sersic profiles. We 
therefore used the Bondi power from the best-fit /3 -model and its 
associated errors. 



The Kendall's t rank correlation was used to determine if 

these two measures of cavity power are significantly correlated with 

the Bondi accretion rate. For our estimates of cavity power, we find 

no significant correlation with T = 0.2. For the estimates of cav- 

I 1 i ~i 

ity power from Allen et al. (2006), we calculate T = 0.7 and reject 

the null hypothesis of no correlation at 95% confidence but not at 

99% confidence. This analysis therefore suggests weaker evidence 

for a trend between the cavity power and Bondi accretion power, 

primarily due to the uncertainty in estimates of the cavity volumes. 



4 DISCUSSION 

X-ray central point sources appear to be common in BCGs host- 
ing X-ray cavities. We find a detection fraction of ~ 50 per cent 
for the BCGs in our sample. The majority of these sources are ra- 
diatively inefficient with required average accretion rates of only 
10 -5 — 10 _2 MEdd- The nuclear X-ray luminosity for these sources 
was observed to correlate with the AGN cavity power, which is 
surprising given the vastly different timescales for these quan- 
tities. Cavity power is averaged over the bubble ages, typically 
10 7 — 10 8 yr, while the nuclear X-ray luminosity is an instantaneous 
measurement and we have shown that this can vary significantly on 
shorter timescales of months to years. The scatter in this correla- 
tion covers over three orders of magnitude. A significant fraction 
of this scatter is likely due to X-ray variability but absorption and 
uncertainty in the cavity power estimates will also contribute. 

The interpretation of these results is complicated by the un- 
certainty in the origin of the nuclear X-ray emission. The X-ray 
emission may originate from the accretion disk corona, from the 
base of a parsec- scale jet or a combination of the two although 
another origin is also possible. However, the nuclear X-ray emis- 
sion is generally considered a probe of accreti on power, whether it 
is from the accretion flow or f rom the jet (eg. iFalcke & Biermarir] 
1 19951 : iHeinz & Sunvaevll2003h . Therefore, the observed L x - P C m 
correlation suggests the accretion power roughly scales with the 
cavity power over long timescales with the large scatter reflecting 
variability on shorter timescales. 



4.1 Duty cycle of activity 

It is also not clear if systems with only upper limits on the point 
source luminosity are simply faint or currently 'off. As shown in 
Fig. |4] whether the upper limits form the faint end of the detected 
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Figure 10. Deprojected temperature and electron density profiles of a subset of the cluster sample for which the cluster properties can be resolved at radii 
within an order of magnitude of the Bondi radius (shown by the vertical dashed line). 
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source population or are instead a different population of 'off sys- 
tems can have a significant impact on the correlation's slope and 
scatter. Approximately ~ 50 per cent of the sample do not have de- 
tected nuclear X-ray emission. It is likely that at least some sources 
are simply a little too faint to be detected, particularly if they are 
embedded in bright cluster emission (see section [3^21 ). These ob- 
jects may therefore still be consistent with the observed P cav — Lx 
correlation. However, objects such as MS0735 and Sersic 159 have 
upper limits on their radiative luminosities a factor of thousand be- 
low that expected from this trend and are effectively 'off. 

Irllavacek-Larrondo & Fabianl d201 lh considered a sample of 
highly radiatively inefficient nuclei in clusters with powerful AGN 
outbursts, including MS0735, and suggested several explanations 
including absorption and variability. We have found significant 
variability for several sources in a subset of our sample, which 
could indicate a cycle of activity, but don't find absorption to be 
as important. Although we caution that this sample is by no means 
complete, the fraction of detections to non-detections indicates a 
duty cycle of at least ~ 50 per cent in systems with recent AGN 
outbursts. If we con sider only the 3 1 sources that overlap with the 
Raff ertvetaD l2006) sample, we find a similar detection fraction of 



at least 40 per cent. This suggests that roughly half of all systems 
undergoing an AGN outburst in the last ~ 1Q 8 yr hav e evidence 
of ongoing accretion. iMendvgraL Jones & Dolad(l2012h found that 
simulations with a jet duty cycle of 50 per cent, cycling on and off 
with a 26 Myr period, produced multiple cavity pairs with a sim- 
ilar morphology to observations ( see also lO'Neill & Jonesll2010l : 
Mendygrai rp Neill & Jonesll201 3). For complete samples of clus- 
ters, the fraction with detected X-ray cavities implies a duty cycle 
of at least ^ 60 — 70 pe r cent ( Dunn & Fabi an 2006; Birz an et al] 
l2009l : lBirzanetalJl2012h . 

4.2 Radiative efficiency and evidence for a transition 
luminosity? 

Although it is not clear if the upper limits form the faint end of 
the detected source population or are a separate population of 'off 
systems, the best fit models for these two possibilities have a con- 
sistent slope in Fig. 01 (left). This slope shows an increase in ra- 
diative efficiency with the mean accretion rate. The quasars in- 
cluded for comparison form an extension of this trend from cav- 
ity power-dominated to radiation-dominated sources. Studies have 
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Figu re 11. Bond i accre tion power versus cavity power. Left: cavity power from this analysis. Right: cavity power from l Allen et al.l ( E006h . The best-fit relation 
from I Allen et al. (2006) is shown as a dashed line. 



also shown that the radio loudness of low luminosity AGN to lu- 
minous quasars is inversely correlated with the mass accretio n rate 
reg. lHoll2002l : lTerashima & WilsorJEoollPanessa et al.ll2007h . Su- 
permassive black holes appear to become more efficient at releas- 
ing energy through jets as their accretion rate drops. Hlavacek- 
Larrondo et al. (submitted) also find strong evolution in the nuclear 
X-ray luminosities of SMBHS hosted by BCGs such that the frac- 
tion of BCGs with radiatively-efficient nuclei is decreasing over 
time. 

Observational evidence suggests that the accretion pro- 
cess is largely similar for both stellar mass and supermassive 
black holes and therefore we could potentially use stud- 
ies of X-ray binaries to u nderst a nd accretion in AGN (eg. 
Macc arone. Gallo & Fended 12003 ; iMerloni. Heinz & di Matted 



minosity of the point source and scaled by the Eddington accretion 
rate, 



20031 : lFalck~K ording & Markoffl 120041) . X-rav binaries are 



broadly classified into low-hard and high-soft states, which relate 
to the accretion disk properties a nd variation in the accretion rat e 
can trigger state transitions (eg. iRemillard & McClinto ck 2006). 
In the low-hard state, the accretion rate is low, the accretion disk 
is optically thin and radiatively inefficient. The mechanical power 
of the radio jet dominates over the radiative power and the X-ray 
and radio fluxes are cor related (eg. iGallo. Fender & Poolevll2003l ; 
Fender & Belloni 2004). Observations of X-ray binaries have 
shown that as the accretion rate rises above ~ 0.01 — OAM^d the 
source makes a spectral transition from the low-hard to the high- 
soft state (eg. lNowakiri995l : bone. Gierlinski & Kubotall2007h . In 
this state the X-ray emission is dominated by an optically thick, 
geometrically thin accretion disk and the radio emis sion drops 
dram atically suggesting the outflow is suppressed feg. lFender et al.1 
1 1 999b . 

Fig. Q2] shows the radiative and cavity power output as a func- 
tion of the required mean accretion rate for our AGN sample, where 
all quantities are scaled by the Eddington rate. The mean accretion 
rate was calculated from the cavity power plus the bolometric lu- 



M _ (Pcav + ^bo 



^Edd 



^Edd 



(4) 



Lbo\ was calculated as shown in section [37TI for the low luminosity 
AGN and was taken from the literature for the quasars. Note that 
for most of the sources considered L\, \ is insignificant compared 
to P cav and the required mean accretion rate is dictated by the cav- 
ity power. The quasars are the obvious exceptions. There are two 
points for each source on the plot showing both the cavity power 
and the radiative power. For sources where the radiative power or 
the cavity power dominates the output, the corresponding points 
will, by definition, lie on a line of equality between Power/L E( jd 
and M/Mndd- This produces a clear line of points along y = x in 

Fig.na 

The illustrative model from Churaz ov et al. 1 d2005h of a change 
from a radiatively inefficient, outflow dominated mode to a ra- 
diation dominated mode has been shown for comparison in Fig. 
[T2l Fig. Q2] shows a trend of increasing radiative efficiency with 
mean accretion rate (see also Fig.|4]). The radiative and mechanical 
power outputs converge and become comparable at an Eddington 
rate of a few per cent. For accretion rates below ~ 0.1M E( jd the 
cavity power dominates over the radiative output, which is a fac- 
tor of 10 — 1000 times lower. Above ~ 0.1M£dd> a transition ap- 
parently occurs where mechanical power drops suddenly and the 
radiative power strongly dominates. This strong transition is seen 
in three objects: H1821+643, IRAS 09 104+4 109 and 3C 186. These 
are quasars in the centres of galaxy clusters, few are known but they 
show this intriguing and potentially very important effect where 
they transition from mechanically dominated to radiation domi- 
nated AGN. This is precisely the behaviour expected when the ac- 
cretion rate increases and an object transitions from an ADAF to a 
geometrically thin and optically thick disk. The AGN in our sam- 
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Figure 12. The required mean accretion rate scaled by the Eddington rate, 
M/MEdd > plotted against the cavity power (blue circles) and the radia- 
tive power (red triangles) scaled by the Eddington luminosity. Therefore, 
there are two points for each source. Detected nuclear X-ray point sources 
are shown by the solid symbols and upper limits are shown by the open 
symbols. Quasar sources are labelled. The radiation and outflow model 
lines are illustrative only and show a transition from outflow domination 
at low accretion rates to radiative domination at high accretion rates (from 
IChurazov et all2005l. Fig. 1). 



pie therefore appear to show the same qualitative behaviour with 
variation in accretion rate found for stellar mass black holes. 

However, there is significant scatter in the trend for radia- 
tively inefficient sources, some of which may be due to variabil- 
ity in the X-ray flux. There are also outliers on the plot, notably 
MS0735 and A2390. A2390 may have overestimated cavity power 
and thus could move to the l eft. MS0735 was note d as anoma- 
lous to Churazov's scenario in llChurazov et al.1 d2005h . In this case 
the i ssue could be related t o powering by the spin of the black 
hole (IMcNamara et alj|2009) where the spin energy is tapped more 
effici e ntly than mc 2 (eg. McNamara, Rohanizadegan & Nulsenl 



20111: iTch ekhovsk ov. Naravan & McKinnevj 1201 ll : ICad 1201 ll : 



McKinney, Tchekhovskoy & Bl andfordf 



2012). This would imply 



greater jet power per accreted mass than objects powered directly 
by accretion, moving it to the left in Fig.[l2l It is also possible that 
the unknown value of £, the conversion efficiency between mass 
and energy P = eMc 2 , is a large source of scatter, particularly as 
it is applied to mechanical power. We have assumed in eq.Blthat £ 
is tied between the radiation and cavity power and divides out but 
this is of course not necessarily true. Nevertheless, the increasing 
nuclear brightness relative to mechanical power is solid. And the 
transition to quasars does depend on power output and by infer- 
ence, M. This picture is also a simplification of stellar mass black 
hole state transitions. Observed transitions from the low-hard to the 
high- soft state in X-ray binaries are accompanied by an intense and 
rapid radi o outburst, which has no obvio us analogy in our AGN 
model (eg. Fender, Belloni & Gallo 200i). 



4.3 Accretion power 

The source of fuel for the observed AGN activity has been 
the subject of considerable debate. There is, in general, suffi- 
cient cold gas in BCGs to fuel the range of observed jet powers 
tedgdl200ll : ISalome & Combesll2003l : ISokerll2008l : iDonahue et ail 
120111) . However, Bondi accretion directly from the cluster's hot 
atmosphere is appealing because it can provide both a steady 
fuel supply and a simple feedback mechanism. Although whilst 
the gas density is high enough to supply sufficie nt fuel to 
low p ower jet systems through Bondi accretion (eg. lAUen et all 
120061 ). this is difficult to achieve for high power jets (> 



e rg s 1 ; Raffertv et al] 20061: lHar dcastle. E vans & Crostonl 



120071 : iMcNamara. Rohan izadegan &Nulsen 201 1|). Using a sam 
pie of nearby systems, I Allen et al.l (2006) found a correlation be- 
tween the cavity power and Bondi accretion rate suggesting that a 
few per cent of the rest mass energy of material crossing the Bondi 
radius emerges in the jets. 

In our analysis of 13 systems, including the lAllen et alJ (f20Q6) 
sample, we found a significantly larger scatter in the correlation 
between the cavity power and Bondi accretion rate. This was pri- 
marily due to differences in our calculation of the density at the 
Bondi radius and estimates of the cavity power. Calculation of the 
density at the Bondi radius required an extrapolation of this profile 
over an order of magnitude in radius for each object in our sample. 
We used three different model density profiles for the extrapolation 
and found each provided a similarly good fit, with the exception 
of the Centaurus cluster where a /3 -model was significantly better. 
This fit was extended to cover a few kpc rather than just the central 
points as cavity substructure close to the centre affected the depro- 
jection of the density profile. The use of three equally plausible 
extrapolations of the density profile produced different estimates 
of the density at the Bondi radi us and increased th e scatter in the 
Bondi accretion rate compared to lAllen et alJ d2006h . Deeper obser- 
vations of several clusters, showing the cavity extent more clearly, 
significantly altered the cavity power measured for those objects. 
We have therefore found weaker evidence for a trend between the 
cavity power and Bondi accretion rate. 

Bondi accretion is energetically a plausible mech- 
anism for fuelling the lower-powered radio sources in 
our sample. However, it is insufficient to fuel the most 
powerful systems fHardcastle. Evans & Crostonl [2007; 
IMcNamara. Rohanizadegan & Nulsenl l201lh . There are also 
theoretical issues with Bondi accretion that include the ability to 
shed angular momentum and the zero central pressure requirement 
(Pro ga & Begelmanl 12003: IPizzolato & Sokerl I2005L l2010l : ISokerl 
120081: lNaravan & Fabian 201 lb . For a more complete discussion 
see IMcNamara & Nulsenl d2012h . 

In view of the high column densities found for many objects 
in this sample, which are consistent with significant levels of cold 
circumnuclear gas, and the prevalence of cold mol ecular gas in cD 
galaxies (eg. Edge 2001; Sal ome & Combes 120031) we suggest that 
cold gas fuelling is a likely source of accretion power in these ob- 
jects. Nevertheless we cannot rule out or exclude Bondi accretion, 
which could play a significant role, particularly in low power jets 
dAUen et alJl2006h . 



4.4 Nuclear X-ray emission mechanism 

The origin of the observed nuclear X-ray emission is not cur- 
rently understood. Observed correlations between the X-ray and 
radio core luminosities provide the strongest support for a 
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non-t h ermal jet-rela t ed or i gin (eg. iFab biano, Gioia & Trin chieril 
Il989l : ICanosa et alJ 1 19991 : iHardcastle & Worralll Il999h . Radio 
and optical luminosity correlations for FR I nucl e i also sup- 
port this conclusion dChiaberge. Capetti & Celottil [l999) and 
multi-wavelength spectral energy distributions for these sources 
can be modelled by sync hrotron and synchro t ron self-Compton 
emission from a jet ( eg. ICapetti et al. 2002; Yuanetal. 2002; 
Chiaberg eltal] l2003ah . However. 
d2004h show that a significant fraction of sources with strong 
optical jet emission do not have an X-ray component po- 
tentially indicating different physical origins. The detection 
of broad FeKa lines and rapid variability on ks timescales 
favou rs an accretion flow origin (eg. iGliozzi. Sambruna & Brandt! 
l2003h . Radio and X-ray correlations do not necessarily im- 
ply a common origin for the emission as accretion pro- 
cesses and jets are likely to be correlated phenomena (eg. 
Bege lman. Blandford & Reesl 1 19841). iM erloni, He inz & di Matted 
(20031) and balcke. K ording & Markoffl ( 120041) argued that these 
correlations are part of a 'fundamental plane' linking radio and X- 
ray emission to black hole mass but diff er on whether this reveals 
trends in accretion or jet physics (eg. iKording. Falcke & Corbell 
2006; Hardcastle, Evans & Cros ton 2009; Plotkin et al.ll2012l) . 



4.4.1 Jet origin? 

IWu. Yuan & Caol d2007h analysed the spectral energy distributions 
of eight FR I sources including two of the variable sources in our 
sample, A2052 and M84 (3C 272.1). They found that the emission 
in M84 is dominated by a jet and the ADAF model predicts too hard 
a spectrum at X-ray energies. A2052 appears to have a comparable 
contribution from the jet and the ADAF. These two sources have 



count rate in the Chandra observations of our sample is generally 
not large enough to search the light curve of each individual obser- 
vation for flux variation. We also do not find any significant varia- 
tion in the spectral properties of the three variable sources identified 
in our subsample. 



very different accretion rates: for A2052 M/M Edd = 9+2 x 10 4 
and for M84 M/M Edd = 5±J x 10~ 6 . Yet they both experience sig- 
nificant variations in nuclear flux on timescales of months to years. 

Strong flux variability on timescales of 1 — 2 months is 
seen from both the core and the jet knot HST -1 in M87 (eg. 
lHarris. Biretta & Junorl Il997l : lHarris et alJ l2006h . From 2000 to 
2009, HST-1 is the site of a massive X-ray, UV and radio flare. Dur- 
ing this period, the X-ray emission from HST-1 dominat es over the 
nucle us and rises and falls by an order of magni tude (lHarris et al. 
2009). The nuclear variability is characterised by Harr is et aI1l2009l 
as 'flickering' with changes in flux of order a factor of a few over 
timescales of months to years. It is not clear if the nuclear X-ray 
emission is due to the inner unresolved jet or the accretion flow. 
The magnitude and timescales of the X-ray flux variability found 
in A2052, Hydra A and M84 are therefore consistent with that ob- 
served in M87. Interestingly, HST UV observations of the nucleus 
in A2052 find the luminosity in creased by a factor of ten from 1994 
to 1999 dChiaberge et al.l l2002). This period was then followed by 
a decrease in the X-ray luminosity by a similar factor from 2000 
to 2010, which could indicate a flaring event similar to that experi- 
ence d by HST-1. 

Sambruna et al 1 (120031) found variability with XMM-Newton 
on 3 — 5 ks timescales in the FR I radio galaxy NGC4261, 
which is also part of our sample. For an ADAF, the X-ray emis- 
sion is radiated from a relatively large volume an d variation 
is expected on timescales l onger than around a day (Pta k et all 
Il998l ; Flbrashima et al.ll2002l) . The observed variability timescale 
in NGC4261 is around two orders of magnitude shorter than the 
ADAF light crossing time suggesting that the variable component 
is more likely to be associated with the inner jet. Unfortunately, the 



ADAF models predict trends between the nuclear radio and X- 
ray lu minosities that scale as LrocL^ 6 and MocL^ 5 jYi & Bo ughn 
1998). Assuming MocP cav , this scaling is consistent with the slope 



shown in Fig. [4] Although the slope between Lr and Lx shown in 
Fig. [9] appears to be steeper than LrocL^ 6 , the lack of a flux-flux 
correlation and the large uncertainties suggest this trend is unre- 
liable (see section [3T6j>. It is therefore worthwhile to consider the 
consequences of emission from an ADAF. 

Fig. [13] compares the X-ray point source luminosities with 
predictions from ADAF models at different accretion rates 
(Merloni, Heinz & di Matteo 2003). The nuclear X-ray luminosity 
scales very close to linearly with black hole mass allowing us to 
scale up these models to the required ~ 10 9 M . For low accretion 
rates, the 2 — 10 keV emission includes inverse Compton scattering 
of soft synchrotron or disk photons and a bremsstrahlung compo- 
nent at higher energies. At higher accretion rates, the Compton- 
scattered component dominates as the optical depth rises and cool- 
ing processes become more efficient. The exact scaling of L2-10 keV 
with M/MEdd W1 ll depend on the parameters chosen for the model, 
such as the visc osity, magnetic pressure and electro n heating frac- 
tion. Therefore jMerloni. Heinz & di Matteol (f2003) obtain a single 
powerlaw fit ^2-lOkeV ^ (M / ^M Edd ) 2,3 from their ADAF models 
for comparison with observational data. 

The required mean accretion rate was calculated from the cav- 
ity power plus the radiative power as shown in eq.|4] For most of the 
sources considered L^ \ is insignificant compared to P cav and the 
required mean accretion rate is dictated by the cavity power. The 
three quasars are the obvious exceptions to this. The majority of the 
observed sources are consistent with the emission expected from 
an ADAF for black hole masses from 5 x 1O 8 M to 5 x 1O 9 M . 
Four sources in our sample have more accurate dynamical black 
hole masses and therefore provide a more reliable constraint when 
compared with ADAF model predictions. Cygnus A and NGC4261 
have mass accretion rates a factor of 3 — 5 smaller than the ADAF 
model predictions for their respective black hole masses. The dif- 
ference is even greater for M84 and M87 with over an order of mag- 
nitude and close to two orders of magnitude discrepancy, respec- 
tively. There are several possible reasons for this. The cavity power 
in M84 is particularly difficult to estimate as the outburst appears to 
have blown out most of the X-ray atmosphere. The cavity volume 
and surrounding pressure may therefore have been underestimated. 
The total mechanical power in M87 has been significantly under- 
estimated as the shock produces an additional 2.4 x 10 43 erg s _1 , 
which is four times greater than the cavity power. This would cause 
M84 and M87 to move to the right in Fig. [13] and closer to the 
ADAF models. So this translates to a similar increase in the mean 
accretion rate. The nuclear X-ray emission from M84 and M87 is 
therefore likely to be consistent with an ADAF given the large un- 
certainties in the ADAF models. As previously discussed in section 
12.41 the cavity power is likely to have been underestimated for the 
majority of the systems in this sample and this will tend to move 
points to the right in Fig. [13] 

We therefore conclude that it is plausible that the X-ray point 
source emission is due to an ADAF but we cannot distinguish be- 
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Figure 13. Point source X-ray luminosity versus the inferred accretion rate 
scaled by the Eddington rate . The lines are ADAF model predictions from 
iMerloni. Heinz & di Matted 120031) for black hole masses of 5 x 1O 8 M 
(dashed), 1 x 1O 9 M (solid) and 5 x 1O 9 M (dash-dotted). The variable 
point sources are shown by the blue stars. Sources with dynamical black 
hole masses are shown by the open symbols (Cygnus A, M84, M87 and 
NGC4261). 



tween this and a jet origin with the available data. Given the lack 
of a clear trend between the nuclear radio and X-ray flux it is likely 
that further progress on this problem will require modelling o f the 
AGN spectral energy distribution feg. lWu. Yuan & Cao 2007). 



5 CONCLUSIONS 

Using archival Chandra observations of 57 BCGs, we have inves- 
tigated the relationship between nuclear X-ray emission and AGN 
radio jet (cavity) power. Although this is not a complete sample 
of objects, we find that nuclear X-ray emission is common with 
roughly half of the sample hosting a detectable X-ray point source. 
Assuming nuclear X-ray emission indicates active accretion, our 
study implies that the AGN in systems with recent outbursts are 
'on' at least 50 per cent of the time. Furthermore, we examine the 
correlation between the nuclear X-ray luminosity and the average 
accretion rate determined from the energy required to inflate the X- 
ray cavities. This correlation is consistent with the hypothesis that 
the nuclear X-ray emission traces active accretion in these systems. 
The majority of the sources in this sample are radiatively inefficient 
with required mean accretion rates of only 10 -5 — 10 _2 MEdd- The 
nuclear X-ray sources become more luminous compared to the cav- 
ity power as the average accretion rate increases. The nuclear X-ray 
emission exceeds the cavity power when the average accretion rate 
rises above a few percent of the Eddington rate, where the AGN 
power output appears to transition from cavity power-dominated to 
radiation-dominated in the three BCGs hosting quasars. 

A small subset of the clusters in our sample had multiple 
archival Chandra observations of sufficient depth to search for 
variability in the nuclear X-ray flux. We found that A2052, M84 
and Hydra A were significantly varying by factors of 2 — 10 on 



timescales of 6 months to ten years. Despite the large variations in 
flux, we did not find significant change in the shape of the nuclear 
spectra. This analysis is generally limited by the availability of suit- 
ably spaced observations of sufficient depth in the Chandra archive 
but suggests that a significant fraction of AGN in BCGs may be 
varying on timescales of months to a few years. This variability is 
likely to be a significant source of the large scatter in the observed 
correlation between the nuclear luminosity and cavity power. Our 
results suggest that the accretion power roughly scales with the cav- 
ity power over long timescales with the large scatter reflecting the 
variability on shorter timescales. 

The interpretation of these results is complicated by the un- 
certainty of the nuclear X-ray emission origin. This emission may 
originate from the accretion disk corona, from the base of a jet or a 
combination of both, although other mechanisms are also possible. 
We discuss the similarity in magnitude and timescale of the X-ray 
variability found in A2052, M84 and Hydra A to that observed from 
both the core and jet knot HST-1 in M87. We also show that the 
X-ray nuclear luminosity and required mean accretion rate of the 
systems analysed are consistent with the predictions from ADAF 
models. We conclude that an ADAF is a plausible origin of the X- 
ray point source emission but we cannot distinguish between this 
and a jet origin with the available data. 

We have also considered the longstanding problem of whether 
jets are powered by the accretion of cold circumnuclear gas or 
accretion from the hot keV atmosphere. For a subsample of 13 
nearby systems, the Bondi accretion rate was calculated using three 
equally plausible model extrapolations of the cluster density profile 
to the Bondi radius. The results suggest weaker evidence for a trend 
between the cavity power and the Bondi accretion rate, primarily 
due to the uncertainty in the cavity volumes. Cold gas fuelling may 
therefore be a more likely source of accretion power given the high 
column densities found for many objects in our sample, which are 
consistent with significant quantities of cold circumnuclear gas, and 
the prevalance of cold molecular gas in BCGs. However, we cannot 
rule out Bondi accretion, which may play a significant role, partic- 
ularly in low power jets. 
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Table 1. Sample properties and point source fluxes. 



Target Redshift Obs. ID Aimpoint Exposure 42 nu b «h,z T ^p,2-i0keV c ^S,2-I0kev c References 











(ks) 


(10 22 cm- 2 ) 


(10 22 cm" 2 ) 




(10- 14 erg cm -2 s _1 ) 


(10- 14 erg cm -2 s _1 ) 




2A0335+096 


0.0349 


7939 


S3 


49.5 


0.2218 





1.9 


< 2.3 


< 0.1 


[1] 


3C295 


0.4641 


2254 


13 


75.6 


0.0134 


43+!? 


1.9 


76 ±4 


55 ±3 


[2], [3] 


3C388 


0.0917 


5295 


13 


26.3 


0.0555 





1.9 


4.0 ±0.3 


2± 1 


[4],[1] 


3C401 


0.2011 


4370 


S3 


22.7 


0.0582 





1.9 


3.7±0.3 


3.7±0.3 


[2],[1] 


4C55.16 


0.2412 


4940 


S3 


64.5 


0.0429 





1 ^^+0 09 


7.3±0.3 




[5],[1] 


A85 


0.0551 


904 


10 


38.2 


0.0278 





1.9 


< 2.2 


< 2.1 


[1] 


A133 


0.0566 


9897 


13 


67.9 


0.0159 





1.9 


< 0.8 


< 0.2 


[1] 


A262 


0.0166 


7921 


S3 


108.6 


0.0567 





1.9 


< 0.49 


< 0.04 


[1] 


A AHQ 
A4/0 


U.Uool 


1 C*C*C\ 

loo9 




39.3 


U.z/ /o 


rv 

u 


1.9 


^ o o 

< 3.2 


< 1.4 


r cs\ r 1 1 

[6],[1] 


A611 


0.2880 


3194 


S3 


32.0 


0.0447 





1.9 


3.0±0.2 




~ 


A1795 


0.0625 


10900 


S3 


15.8 


0.0119 





1.9 


< 1.8 


< 1.6 


[1] 


A1835 


0.2523 


6880 


13 


109.2 


0.0204 





1.9 


< 3.1 


< 1.2 


[6],[1] 


A2029 


0.0773 


4977 


S3 


74.7 


0.0325 





1.9 


< 4.5 


< 1.4 


[1] 


A2052 


0.0351 


5807 


S3 


123.8 


0.0272 





1.9 


3.8±0.1 


a i -l-u.z 


[7],[1] 


A2199 


0.0302 


10748 


13 


40.6 


0.0089 





1.9 


3.1±0.3 


0.9±0.6 


[8],[1] 


A2390 


0.2280 


4193 


S3 


70.6 


0.0768* 





1.9 


1.4±0.1 


A ' u -0.2 


- 


A2597 


0.0852 


7329 


S3 


56.3 


0.0248 





1.9 


< 2.3 


< 1.3 


[1] 


A2667 


0.2300 


2214 


S3 


8.4 


0.0173 





1.9 


2.4 ±0.4 


2.8±0.4 


~ 


A4059 


0.0475 


5785 


S3 


85.9 


0.0121 





1.9 


< 0.72 


< 0.08 


[1] 


Centaurus 


0.0114 


4954 


S3 


87.2 


0.0948* 





1.9 


< 1.6 


< 0.10 


[6],[1] 


Cygnus A 


0.0561 


1707 


S3 


9.2 


0.272 


19.6±0.7 


1.9 


3160±70 


2560 ± 60 


[6],[1] 


HCG62 


0.0137 


10462 


S3 


65.1 


0.0332 





1.9 


< 0.67 


< 0.03 


[1] 


Hercules A 


0.1540 


6257 


S3 


48.5 


0.0618 





1.9 


< 0.94 


< 0.53 


[1] 


Hydra A 


0.0549 


4970 


S3 


94.6 


0.0468 


3.5±0.3 


1.9 


24 ±1 


. A ,-_|_o 8 


[9],[1] 


M84 


0.0035 


5908 


S3 


45.1 


0.0299 


~ , _ _i_n rn 
u,iZ -0.02 


2.13±0.09 


17.5±0.4 


15 -3 


[10],[1] 


M87 


0.0044 


1808 


S3 


12.8 


0.0194 


0.06 ±0.01 


2.37 ±0.07 


67±1 


66±4 


[H],[l] 


M89 


0.0011 


2072 


S3 


53.4 


0.0262 


0.092 ±0.03 


2.3±0.1 


5.5±0.2 


3±2 


[12],[13] 


MKW3S 


0.0450 


900 


13 


51.9 


0.0268 





1.9 


< 1.0 


< 0.2 


[1] 


A/TQn7^^ £_i_7ZL91 


U.ZIOU 




Lj 


1 JJ.O 


U.UjZo 


n 
u 


1 Q 


/HQS 


<. U. 1Z 


L14J 


NGC507 


0.0165 


2882 


13 


40.0 


0.0525 





1.9 


< 1.1 


<0.1 


[15] 


NGC1316 


0.0059 


2022 


S3 


20.2 


0.0240 





1.9 


<3.0 


<2.2 


[15] 


NGC1600 


0.0156 


4283 


S3 


20.5 


0.0349 





1.9 


<0.6 


<0.2 


[15] 


NGC4261 


0.0075 


9569 


S3 


99.9 


0.0175 


3 +0 - 6 


0.9±0.2 


126 ±2 


74±2 


[16],[17] 


NGC4472 


0.0033 


11274 


S3 


39.7 


0.0153 





1.9 


<2.4 


<0.1 


[15] 


NGC4636 


0.0031 


3926 


13 


67.8 


0.0190 





1.9 


<0.9 


<0.3 


[15] 


NGC4782 


0.0154 


3220 


S3 


49.3 


0.0337 





1.9 


0.84 ±0.09 


0.7±0.1 


[18],[15] 


NGC5044 


0.0093 


9399 


S3 


82.5 


0.0487 


n 9 o+0.08 


1.9 


2.6±0.2 


1 4+0.4 


[19],[15] 



Table 1 (cont'd) 



Target 


Redshift 


Obs. ID 


Aimpoint 


Exposure* 3 
(ks) 


(10 22 cm -2 ) 


At +L,Z 

(10 22 cm -2 ) 


r 


J7 C 

1 r,z— 1U keV 

(10~ 14 erg cm -2 s _1 ) 


J7 C 

1 5,2— 1U kev 

(10~ 14 erg cm -2 s _1 ) 


References 


NGC5813 


0.0066 


9517 


S3 


98.8 


0.0918* 





1.9 


< 1.02 


<0.3 


[20] 


NGC5846 


0.0057 


7923 


13 


84.7 


0.0429 





1.9 


<0.72 


<0.05 


[15] 


NGC6269 


0.0348 


4972 


13 


35.3 


0.0510 





1.9 


<2.4 


<0.08 


[15] 


NGC6338 


0.0274 


4194 


13 


44.5 


0.0223 





1.9 


<2.2 


< 1.4 


[15] 


PKS 0745-191 


0.1028 


12881 


S3 


116.0 


0.38 





1.9 


4.8±0.3 


2 6 +0 - 2 


[1] 


PKS 1404-267 


0.0230 


12884 


S3 


83.5 


0.0578* 


0.04 ±0.02 


2.01 ±0.06 


21.4±0.4 


19.8±0.7 


[21],[1] 


RBS797 


0.3540 


7902 


S3 


38.3 


0.0228 




1.4±0.2 


61±2 


41 ±2 


[22] 


RXCJ0352.9+1941 


0.1090 


10466 


S3 


27.2 


0.1440* 


3.6±0.4 


1.9 


36±2 


32±2 




RXCJ1459.4-1811 


0.2357 


9428 


S3 


39.5 


0.0736 


c 4+i. u 
J -^-0.8 


1.9 


21 ±1 


13±1 




RXCJ1524.2-3154 


0.1028 


9401 


S3 


40.9 


0.1107* 


0.5 ±0.2 


1.9 


5.7±0.5 


3.2±0.4 




Ix/VV^J 1 JJO.j-14-lu 


O 0Q70 


QZL09 


OJ 


^(\ o 

jD.U 


o io^o 


ic+5 


1 Q 


H-l ZC J 


99-1-9 
zz m z 




Sersic 159-03 


0.0580 


11758 


13 


93.3 


0.0114 





1.9 


< 1.0 


<0.07 


[1] 


UGC408 


0.0147 


11389 


S3 


93.8 


0.0256 


O oa+0.06 
U - ZJ _0.05 


2.2±0.2 


3.3±0.1 


2.8±0.2 


[20] 


Zw348 


0.2535 


10465 


S3 


48.9 


0.0250 




1.9 


4.4 ±0.4 


2.9±0.4 




Zw2089 


0.2350 


10463 


S3 


38.6 


0.0286 


22±2 


1.9 


203 ±7 


127+^ 




Zw2701 


0.2150 


12903 


S3 


95.8 


0.0075 





1.9 


<0.53 


<0.3 


[6],[23] 


Zw3146 


0.2906 


9371 


13 


34.5 


0.0246 





1.9 


<5.0 


<3.3 


[6],[1] 



Note. — ^Final cleaned exposure times. ^Column densities marked with * were found to be higher than the Galac tic value measured b y iKalberla et alJ d2005h (see section [231 . 
^E rrors and upper limits on the phot omet ric Fp ^-lokey and spec troscopic Fs T 2-iokeV flux values are l a. References: [1] Raffer tv et alJ (120061) . [2]|Hardcastle, Evans & Croston ( 2006), 
[3] IHlavacek-Larrondo et al] (120121). T41 lEvans et all (120061) . f51 IHlavacek-Larrondo et"aD (1201 ll) \6^ IHlavacek-Larrondo & Fabianl d201lh. ?71 iBlanton Sarazin & McNamarJ (12003), 
[8] iDi Matteo et all (1200 ll) . T91 iMcNamara et al.1 (l200ol). Q0 1 iBalmaverde. Capetti & Grandil ( l2006h. fill iDi Matteo et al.1 j2003l), [12] IMachacek et all d2QQ6b, [13] I Allen et all (12 006). 
[14] McNamara et al. i n pre p Q51 ICavagnolo et al I feOlCh. Q61 IChiaberge et al I (l2003bh . Q71 lO' Sullivan etaD feOllbh . 
lO'Sullivan etaD (1201 lal) . [21] Johnst one et al.1 (l2005h . [221 ICavagnolo et al.N201ll) . [23] McNamara et al. in prep 



[18] IMachacek et all (2007), [19] iDavid et all d2QQ9h . [20] 



Oq 



Oq 
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Table 2. Bondi parameters for the selected subsample of targets. 



Target D L T(r B ) r B n e (r B ) M B Pq P cav 

(Mpc) (keV) (kpc) (cm" 3 ) (Movr" 1 ) (lO^ergs" 1 ) (lO^ergs" 1 ) 



A2199 


132.3 


1.10±0.07 


0.020 ±0.001 


H+o.s 


006+ 0004 
u.uuo_ 0003 


4±2 


9tt 


Centaurus 


49.3 


0.41 ±0.04 


0.028 ±0.002 


0.4+U! 

^ -0.2 


0.002+^02 

^ ^ -0.001 


0.51 ±0.06 


14+] 

-4 


HCG62 


59.3 


0.31 ±0.03 


0.058 ±0.006 


3+ 03 


006+0- 006 
u.uuo_ 0000 


5±2 


6±j 


M84 
M87 


17.0 
17.0 


0.34 ±0.01 
0.52 ±0.04 


0.084 ±0.003 
0.42 ±0.04 


041 +o.io 
u - 4i -o.oi 
0.23 ±0.05 


020+ 005 
u - uzu -0.ooi 

0.35 ±0.07 


13±1 
200 ±40 




M89 
NGC507 


17.0 
71.6 


0.35 ±0.01 
0.49 ±0.02 


0.041 ±0.002 
0.048 ±0.002 


0.8l j 


007+0- 008 
u.uu/_ 0006 

02+0- 04 


5±4 
20 ±10 


19+' 4 


NGC1316 
NGC4472 


25.4 
17.0 


0.343 ±0.007 
0.372 ±0.007 


0.0267 ±0.0005 
0.061 ±0.001 


i + 1.6 
i - 3 -0.9 

0.5 ±0.2 


006+ 008 
u.uuo_ 0004 

nm 4 +0.006 


5±4 
9±3 


7+°- 5 


NGC4636 
NGC5044 
NGC5813 
NGC5846 


17.0 
40.1 
28.4 
24.5 


0.23 ±0.02 
0.31 ±0.01 
0.33 ±0.01 
0.378 ±0.009 


0.028 ±0.002 
0.0290 ±0.0009 
0.0288 ±0.0009 
0.0256 ±0.0006 


20+ 003 
3+°- 5 

0.28_ 009 
4 +0 - 4 


0009+0 001 

u - uuuy -0.0002 

002+ 003 
0.0016^-0010 
0.002^0-002 


0.5±0.1 

2±1 
1.0±0.4 
1.3±1.0 


27+0- 15 
UmA ' -0.09 
i o+0.8 
i - 3 -0.4 

7+°- 6 

1 6+ 12 



Table 3. Point source fluxes and key parameters for each source tested for variability. The photometric point source flux, Fp, and the spectroscopic point 
source flux, Fs, are both given in the 2—10 keV energy band. A spectroscopic flux measurement could not be produced for M84 obs. ID 401 because the 1 ks 
exposure was too short. 



Target 


Obs. ID 


Date 


Aimpoint 


Exposure 
(ks) 


«H,z 

(10 22 cm" 2 ) 


r 


(10" 14 ergcm- 2 s" 1 ) 


(10- 14 ergcm" 2 s- 1 ) 


A2052 


890 


03/09/2000 


S3 


30.5 


< 0.014 


2 06+0- 08 

z - uo -0.06 


12.6 ±0.4 


10±1 




5807 


24/03/2006 


S3 


123.8 


< 0.019 


9 nO+0- 09 


3.8±0.1 


3 l+°- 2 




10879 


05/04/2009 


S3 


80.0 


16+0- 08 

U - AO -0.07 


3.2±0.4 


2.2±0.1 


1.1 ±0.3 




10478 


25/05/2009 


S3 


119.0 


03+0- 05 


2.2±0.2 


2.1 ±0.1 


1.4±0.3 




10477 
10479 


05/06/2009 
09/06/2009 


S3 
S3 


59.0 
63.9 


o.ii±2:S 

0.08^0-08 


9 8+0- 5 

2.8 ±0.4 


2.4±0.1 
2.2±0.1 


1 4+0-3 
1 o+0.3 


A2390 


500 


08/10/2000 


S3 


8.8 





1.9 


0.9 ±0.4 


< 1.1 




4193 


11/09/2003 


S3 


70.6 


l+o- 2 
u - A -o.i 


1.9 


1.4±0.1 


9+°- 1 


Hydra A 


576 


02/11/1999 


S3 


17.4 


4 ^+0-8 


1.9 


36±3 


26±2 




4969 


13/01/2004 


S3 


62.1 


o +0.8 


1 7+°- 4 


27±1 


18±1 




4970 


22/10/2004 


S3 


94.6 


9 9+0.6 
Z ' Z -0.5 


1.2±0.3 


24±1 


14.5^-8 


M84 


401 


20/04/2000 


S3 


1.7 





1.9 


7±1 






803 


19/05/2000 


S3 


25.5 


0.17 ±0.04 


2.1 ±0.1 


8.6 ±0.4 


8.2+? : 4 




5908 


01/05/2005 


S3 


45.1 


12+0- 03 
u>i -0.02 


2.13±0.09 


17.5±0.4 


15* 




6131 


07/11/2005 


S3 


35.8 


17+008 
U>1 7 -0.07 


2.2 ±0.3 


4.9 ±0.3 


4 9+0-3 


NGC5044 


798 


19/03/2000 


S3 


14.3 





1.9 


2.4 ±0.4 


<0.2 




9399 


07/03/2008 


S3 


82.5 


n 90+0.O8 

u - z - 3 -0.07 


1.9 


2.6 ±0.2 




PKS0745 


2427 


16/06/2001 


S3 


17.9 





1.9 


3.0±0.6 


2-ol?:I 




12881 


27/01/2011 


S3 


116.0 


<0.08 


1.8±0.2 


4.8±0.3 


2 6+°- 2 
Z -°-0.4 


PKS1404 


1650 


07/06/2001 


S3 


7.1 


0.07 ±0.05 


2.4 ±0.2 


24±1 


22±2 




12884 


03/01/2011 


S3 


83.5 


0.04 ±0.02 


2.01 ±0.06 


21.4±0.4 


19.8±0.7 
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